Spectral function of the one-dimensional Hubbard model away from half filling 
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We calculate the photoemission spectral function of the one- dimensional Hubbard model away 
from half filling using the dynamical density matrix renormalization group method. An approach 
for calculating momentum-dependent quantities in finite open chains is presented. Comparison with 
exact Bethe Ansatz results demonstrates the unprecedented accuracy of our method. Our results 
show that the photoemission spectrum of the quasi-one-dimensional conductor TTF-TCNQ provides 
evidence for spin-charge separation on the scale of the conduction band width. 
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The Luttinger liquid theory describes the ground state 
and asymptotic low energy properties of one-dimensional 
correlated metals Q. Two characteristics of a Luttinger 
liquid are the absence of quasi-particles predicted by the 
Fermi liquid theory of normal metals and the occur- 
rence of independent spin and charge excitations. In 
principle, these features can be observed in the spec- 
tral function [3, which corresponds to the spectrum 
measured in angle-resolved photoemission spectroscopy 
(ARPES) experiments. In real materials, however, the 
low-energy properties are likely to be governed by three- 
dimensional physics. One-dimensional physics is ob- 
served only above a crossover energy scale, even in the 
most strongly anisotropic materials. Consequently, it has 
proven difficult to observe unambiguous evidence for Lut- 
tinger liquid physics in experiments probing only low- 
energy properties of quasi-one-dimensional conductors. 

A recent ARPES experiment for the quasi-one-dimen- 
sional organic conductor TTF-TCNQ (tetrathiafulvalene 
tetracyanoquinodimethane) has revealed significant dis- 
crepancies from the predictions of Fermi liquid theory 
and conventional electronic structure calculations [j, |5j . 
The experimental spectrum dispersion can be consis- 
tently mapped over the scale of the conduction band 
width onto separated spin and charge excitation bands of 
the one-dimensional Hubbard model away from half 
filling. This is one of the strongest pieces of experimen- 
tal evidence for spin-charge separation and thus for Lut- 
tinger liquid physics in low-dimensional materials. How- 
ever, a direct comparison of the experimental ARPES 
spectrum with the Hubbard model spectral function has 
not been possible yet. 

The Hubbard model was solved exactly 36 years ago 
and the dispersion of its excitation bands can be com- 
puted 0. Nevertheless, the photoemission spectral func- 
tion can only be calculated exactly in the limiting cases of 
noninteracting electrons or infinitely strong electron in- 
teraction llflj and in the low-energy limit described by 
the Luttinger liquid theory. Various numerical methods 



have provided a qualitative picture of spectral functions 
in the Hubbard model but exact diagonalizations 0, 0| 
are limited to too small systems to investigate the ther- 
modynamic limit while other approaches [l^ Il3l| are 
based on various approximations of uncertain accuracy. 

We have determined the photoemission spectral func- 
tion of the one-dimensional Hubbard model with pa- 
rameters appropriate for TTF-TCNQ using the dynam- 
ical density- matrix renormalization group (DDMRG) 
method A novel approach is used to calculate 

momentum-dependent quantities in finite open chains. 
This allows us to investigate large systems almost ex- 
actly and to make a direct comparison of the Hubbard 
model spectral weight distribution with the experimen- 
tal TTF-TCNQ spectrum. To demonstrate the accuracy 
of our method and to identify excitations contributing 
to the photoemission spectral function we compare our 
numerical results with exact Betlic Ansatz results. 

A minimal model to describe the electronic properties 
of TTF-TCNQ is the onc-dimensional Hubbard model 
defined by the Hamiltonian 
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Here c^^, ci^„ are creation and annihilation operators 
for electrons with spin a at site / = 1, . . . , L 

(representing a 7r-type Wannier orbital centered on a 
TCNQ molecule), hi^a = cl^ci^a, and fii = fii^^ -t-nj^. 
Appropriate parameters for TTF-TCNQ are an on-site 
Coulomb repulsion U = 4.9t and a hopping integral 
t = OAeV ^1^. (These values are appropriate for the 
TTF-TCNQ surface, which is probed in ARPES exper- 
iments, not for bulk TTF-TCNQ.) Although the filling 
of the TCNQ band is n = 0.59, we use a slightly dif- 
ferent filling n = 0.6 in our simulations to facilitate the 
finite-size-scaling analysis. For a chain with L sites and 
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N = nL electrons the chemical potential /j, is chosen so 
that Eo{N - 1) = Eo{N + 1), where Eo{N ± 1) is the 
ground state energy with iV±l electrons. Thus the Fermi 
energy is e^? = in the thermodynamic limit L ~> oo. 

The photoemission spectral function A{k,uj) is the 
imaginary part of the one-particle Green's function 
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where j-^o) and Eq are the ground state wavefunction 
and energy of the Hamiltonian This function can 
be calculated for finite broadening rj and system sizes L 
using the dynamical DMRG method The spectral 

properties in the thermodynamic limit can be determined 
using a finite-size-scaling analysis |Q| with an appropri- 
ate broadening ri(L). Here we have used rjL = 9t and 
system sizes up to L = 150 sites. DMRG truncation er- 
rors are negligible for all results presented here (up to 
771 = 400 density-matrix eigenstates have been kept per 
block in our calculations.) 

The operators Ck,a are usually defined using Bloch 
states [i.e., the one-electron eigenstates of the Hamil- 
tonian ^ with periodic boundary conditions ci^L,a — 
Ci^CT in the non-interacting limit {U = 0)]: c/c.o- = 
L^^/^ e^^'^'Q^o- with momentum fc = ^-KzjL for in- 
tegers —L/2 < z < L/2. Since DMRG calculations can 
be performed for much larger systems using open bound- 
ary conditions, it is desirable to extend the definition of 
the spectral function A{k,uj) to that case. Therefore, we 
use the eigenstates of the particlc-in-a-box problem [i.e., 
the one-electron eigenstates of the Hamiltonian on an 
open chain for ?7 = 0] to define the operators 



L 



— sin(fcZ)Q,, 



(3) 



with (quasi-)momentum k = ttz/ {L -|- 1) for integers 1 < 
z < L. Both definitions of Ck,a should be equivalent 
in the thermodynamic limit L — > cxd. Tests for finite 
systems up to L = 32 sites show that both approaches 
are consistent except in the asymptotic Luttinger liquid 
regime [i.e., at low energy (jwj ~ close to the Fermi 
vector kp = TTn/2 {\k — kp\ ^ 1/-^)]- Therefore, open 
chains and the definition © can be used to investigate 
the spectral function A{k,uj). In this paper we present 
only results obtained using this approach. 

Figures n and 121 show the spectral function calculated 
with DDMRG in a chain with i = 90 sites. Since the 
spectrum is symmetric, A{—k,Lu) = A{k,uj), we show 
results for fc > only. Three dispersing features are 
clearly visible in the spectrum for |fc| < kp in Fig^ At 
small binding energy — w there are intense peaks with 
a narrow dispersion (from u; « at fc = ±kp to lo ~ 
—0.5t at fc = 0). This feature corresponds to the spinon 
branches in the Luttinger liquid regime. Note that both 
spinon branches (for fc < and fc > 0) join at fc = and 
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FIG. 1: Line shapes (left) and gray-scale plot (right) of the 
spectral function A{k,uj) for < fc < calculated with a 
broadening 77 = O.lt using DDMRG. 



thus form just one spinon band. At energies w lower than 
the spinon band there is a second spectral feature made 
of peaks with less spectral weight and a wider dispersion 
(from ti> « at fc = ±kp to w w — 1.5t at fc = 0). It 
merges with the spinon band for fc — » ±kp because of 
the finite broadening. This feature corresponds to the 
two holon branches of the Luttinger liquid theory. The 
third spectral feature is made of weaker peaks and has an 
(apparently) inverted dispersion (starting at uj w —1.5t 
for fc = and reaching w ~ —2.2t at fc = ±kp). These 
so-called shadow bands 0] are actually the continuation 
of the holon bands. Thus the second and third features 
correspond to two holon/shadow bands crossing at fc = 0. 
While the spectral weight of the structure associated with 
the spinon and holon bands remains relatively constant 
for all jfcj < kp, the shadow bands rapidly loose intensity 
with increasing fc. 

For |fc| > kp the spectral weight is much lower than 
for |fc| < kp (see Fig. |21). Nevertheless, one can observe 
four dispersive structures in the spectral function. First, 
the shadow band continues from fc — ±kp to ±3kp, 
but its energy increases with jfcj and approaches zero for 

fcj — 3kp. Weaker peaks are also visible at higher en- 
ergy oj than the shadow band for kp < \k\ < 2kp. The 
corresponding binding energy — w increases from about 
zero at fc = zLkp to about 1.7t at k — ±2kp, where this 
second feature meets the shadow band and apparently 
disappears. The third dispersing feature corresponds to 
very weak peaks (not visible on the scale of Fig. [21 with 
energies from w « —3.7t at fc = tt to w —2.2t at 

fcj « 2kp. Note that, despite its weakness, this feature 
corresponds to the spectrum maximum for k k, n. The 
last feature is a sharp drop of the spectral weight at low 
energy. It goes from w — 3.25t at |fc| k kp to lo ~ —A.Qt 
a,t k — TT. We interpret this drop as the lower edge of 
the photoemission spectrum. The little spectral weight 
found at lower energy lo is due to the finite broadening 77 
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FIG. 2: Line shapes (left) and gray-scale plot (right) of the 
spectral function A{k,^) for kp < k < tv calculated with a 
broadening 77 = O.lt using DDMRG. 



used in our DDMRG calculations. Note that these third 
and fourth spectral features are not visible for small |fc| 
because they are too close to the broadened and compa- 
rably much stronger peaks belonging to other structures. 

Figure 13 shows the dispersion uj{k) of the various fea- 
tures found in the DDMRG spectrum for the 90-site 
chain. One clearly sees that the shadow bands are just 
the continuation of the holon bands. The dispersions 
uj{k) should naturally correspond to specific excitation 
bands e(fc) of the Hubbard model. To identify these exci- 
tations we have calculated the excitation energies e(fc) for 
the removal of an electron in the Hubbard model on a 90- 
site chain using the Bethe Ansatz solution |^ . In Fig. El 
we show those excitation bands e(fc) which correspond 
to the dispersing features found in the DDMRG spectral 
function. The excellent quantitative agreement between 
the Bethe Ansatz results calculated for periodic bound- 
ary conditions and our numerical data confirms that the 
open chains used in our DDMRG calculations do not af- 
fect the spectral properties significantly. 

Due to the separation of spin and charge dynamics, 
electron-removal excitations with momentum k are made 
of independent spin and charge excitations with momenta 
kg and kc = k ^ kg, respectively. The spinon band be- 
tween —kp and kp is related to excitations with the low- 
est possible binding energy for fee = and | fcs | < kp. This 
defines the spinon dispersion es(fes), which has a width of 
about 0.5t and gives the spectral onset for |fc| < kp. The 
holon/shadow bands going from —kp to 3k p and from 
—3kp to kp correspond to excitations with the lowest 
possible binding energy for |fes| — kp, < |fec| < ikp, 
and kgkc < 0. This defines the holon dispersion edkc) 
with a width of about 2t. It gives the spectral onset for 
2kp < \k\ < 3kp. The peaks found at low binding energy 
for kp < \k\ < 2k p correspond to secondary holon bands 
made of similar excitations as the holon-shadow bands 
but with parallel spin and charge momenta (fegfcc > 0). 
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FIG. 3: Dispersion ijj{k) of the structures observed in the 
DDMRG spectral function; spinon band (square), holon- 
shadow bands (circle), secondary holon bands (diamond), 
lowest "4fcF" -singlet excitations (plus), and lower (open tri- 
angle) and upper (solid triangle) spectrum edges. Lines show 
dispersions e(k) obtained from the Bethe Ansatz solution. 



They give the spectral onset for kp < \k\ < 2kp. For 
3kp < \k\ < TT this onset corresponds to a secondary 
spinon band with kc — iAkp and \ks\ < kp. 

In Fig. 131 a dashed line shows the dispersion of the 
lowest possible excitations made of one spinon and one 
holon [i.e, the minimum of edkc) + es(fes) for a given 
k = kc + kg]. This lower edge of the spinon-holon contin- 
uum is not related to any feature in the DDMRG spec- 
tral function and one finds spectral weight at lower en- 
ergy u. Therefore, the Hubbard model spectral func- 
tion can not be explained with spinon-holon excitations 
only. Actually, the lower edge of the spectrum follows 
the dispersion of the lowest states made of one spinon 
and a single charge excitation called "4fei?" -singlet exci- 
tation in Ref. 8] . Finally, the very weak peaks found for 
~2t > uj > —At and |fe| > 2kp seem to be related to 
the lowest possible "4fci?"-singlet charge excitations with 
kg — ±kp and kckg > 0. 

In Ref. 01 it was shown that the dispersion of the 
TCNQ related peaks in the ARPES spectrum of TTF- 
TCNQ could be mapped onto excitation bands of a one- 
dimensional Hubbard model. Our DDMRG calculations 
show that the Hubbard model also explains qualitatively 
the experimental spectral weight distribution. (A quan- 
titative comparison is not possible because of the strong 
background contribution in the ARPES data.) The 
ARPES spectrum features labeled (a), (b), and (d) in 
Refs. y, I3I perfectly match the (singular) features found 
in the Hubbard model spectral function (the spinon, 
holon, and shadow bands, respectively). This confirms 
that the ARPES spectrum of TTF-TCNQ shows the sig- 
nature of spin-charge separation over the scale of the con- 
duction band width (of the order of leV). In addition, we 
note that the secondary holon bands (for kp < \k\ < 2kp) 
correspond to a poorly understood spectral feature [la- 
beled (c)] which has been attributed to excitations of the 
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FIG. 4: Spectral functions A{k ~ tt/IO = kp/S, u) calculated 
with DDMRG for system sizes L = 30 (dotted), 60 (dashed), 
and 90 (solid). Inset: scaling of the peak maxima for 0.3 < 
Ti/t < 0.06 (30 < L < 150). Sohd lines are fits. 

TTF band in Refs. HQ. Therefore, we think that this 
spectral feature is not related to the TTF band but is nat- 
urally explained by the TCNQ secondary holon bands, at 
least in the range kp < \k\ < 2kp. 

In the Luttinger liquid theory the spectral functions 
A(k,LS) have singularities \bj — e{k)\~°' for energies e(fc) (x 
\k ± fepj given by the spinon and holon linear disper- 
sions 13, |3| . For a system which is invariant under spin 
rotation the exponents a are related to the Luttinger liq- 
uid parameter Kp through as = (4 — — Kp^)/4: on 
the spinon branch and ac = (6 — Kp — Kp^)/8 on the 
holon branch. The parameter Kp can be calculated in 
the one-dimensional Hubbard model and one finds 
Kp « 0.68 for U = 4.9t and n = 0.6, which corresponds 
to exponents as = 0.46 and ac = 0.48. 

In view of the Luttinger liquid theory results it is nat- 
ural to ask whether the broadened peaks found in our 
DDMRG calculations become singularities of the spec- 
tral function in the thermodynamic limit. To answer this 
question and to estimate the exponents a we have per- 
formed a finite-size-scaling analysis The spectral 
function A{k,uj) is calculated for several system sizes L 
with a broadening scaling as r] — 9t/L. Some spectra for 
fc « tt/IO = kp/S are shown in Fig.^ The scahng of the 
peak maxima A,„ax with rj can then be analyzed (see the 
inset of Fig. 0. If A^ax diverges as 77"" (0 < a < 1) for 
?7 — > 0, the spectral function has a singularity with ex- 
ponent a in the thermodynamic limit. A Landau quasi- 
particle corresponds to a Dirac J-function and thus to a 
peak diverging as ri~^. 

Using this scaling analysis we have found that the 
spinon, holon and shadow band peaks become singular- 
ities in the thermodynamic limit. We have not found 
any diverging peak with an exponent larger than 0.86, 
which confirms the absence of Landau quasi-particles. 



For fc = 7r/10 we have found that the spinon, holon, and 
shadow band exponents are a = 0.78, 0.44, and 0.56, re- 
spectively. For fc = 0, we have obtained a — 0.86 for the 
spinon band and a = 0.70 for the holon/shadow band. 
Therefore, the exponents a are momentum-dependent 
and for finite \k ± kpl they are significantly different 
from the Lut ting er liquid predictions for |fc| kp. A 
recent study has also shown that these exponents 
are strongly fc-dependent. It is not possible to deter- 
mine the exponents a in the asymptotic Luttinger liquid 
regime with DDMRG because the finite-size effects are 
not under control in that limit. 

In summary, we have used a novel approach to com- 
pute the photoemission spectral function of the Hubbard 
model on open chains using DDMRG and explained the 
ARPES spectrum of the organic conductor TTF-TCNQ. 
Our method can easily be extended to other dynamical 
response functions and more complicated models. 
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